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ABSTRACT 



Context. The HsO^ molecule probes the chemistry and the ionization rate of dense circumnuclear gas in galaxies. 

Aims. We use the HjO^ molecule to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar 

medium. 

Methods. Using the JCMT, we have observed the - 2^ 364 GHz line of P-H3O* towards the centers of seven active galaxies. 
Results. We have detected P-H3O+ towards IC 342, NGC 253, NGC 1068, NGC 4418, and NGC 6240. Upper limits were obtained 
for IRAS 15250 and Arp 299. We find large H3O+ abundances (WCHjO+VyVCHj); lO"'*) in all detected galaxies apart from in IC 342 
where it is about one order of magnitude lower We note, however, that uncertainties in ^(HsO^) may be significant due to lack of 
definite information on source size and excitation. We furthermore compare the derived ^(HsO^) with A'(HCO^) and find that the 
H3O+ to HCO+ column density ratio is large in NGC 1068 (24), moderate in NGC 4418 and NGC 253 (4-5), slightly less than unity in 
NGC 6240 (0.7) and lowest in IC 342 (0.2-0.6). We compare our results with models of X-ray and photon dominated regions (XDRs 
and PDRs). 

Conclusions. For IC 342 we find that a starburst PDR chemistry can explain the observed H^O* abundance. For the other galaxies, 
the large HsO^ columns are generally consistent with XDR models. In particular for NGC 1068 the elevated A'(H30^)/A'(HC0^) ratio 
suggests a low column density XDR. For NGC 4418 however, large HC3N abundances are inconsistent with the XDR interpretation. 
An alternative possibility is that HsO^ forms through H2O evaporating off dust grains and reacting with HCO^ in warm, dense gas. 
This scenario could also potentially fit the results for NGC 253. Further studies of the excitation and distribution of HjO^ - as well as 
Herschel observations of water abundances - will help to further constrain the models. 

Key words, galaxies: evolution — galaxies: individual: IC 342, NGC 253, NGC 1068, NGC 4418, NGC 6240, IRAS 15250, Arp 299 
— galaxies: starburst — galaxies: active — radio lines: galaxies — ISM: molecules 



1. Introduction 

Molecular line emission is an important tool for probing the 
highly obscured inner regions of starburst galaxies and buried 
AGNs. Line ratios within and between species help determine 
physical conditions and chemistry of the gas, which provide 
essential clues to the type and evolutionary stage of the nu- 
clear activity. Important extragalactic probes of cloud prop- 
erties include molecules such as HCN, HCO+, HNC , CN , 
and HC3N (e.g. [C ostagliola & Aalto 2010; Krips et a ll 12008 



Gracia-Ca rpio et al. 20 08,; .Loenen et aL ,2008i 



2004t iGao & Solomoril ]2004t lAalto et aLl l2002h that trace the 



dense (nl 10 cm" ) star forming phase of the molecular gas. 
HCN, HNC, HCO+ and CN are all species that can be as- 
sociated both with photon dominated regions (PDRs) (e.g. 



iTielens & HollenbachI 1985h in s tarbursts and X-ray dominated 



regio ns (XDRs) (e.g. iMalonev e t al. 1996; Meiierink & Spaani 
12005) surrounding active galactic nuclei (AGN). In contrast, 
HC3N requires shielded dense gas to survive in significant abun- 
dance since it is destroyed by UV a nd particle radiation (e.g . 
reactions with the ions C^ and He ^' xlPrasad & HuntresslflMt 
[Rodrigue z-Franco et a n [19981 e.g.). These ions are expected to 
be abundant in for example XDRs. Thus HC3N line emission 
may identify galaxies where the starburst is in the early, embed- 
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ded, stage of its evolution. There is still, however, substantial di- 
chotomy in the interpretation of the line emission from the above 
molecules and therefore new tracer species to combine with ex- 
isting information are important. 

The hydronium ion HsO^ is a key species in the oxygen 
chemistry of dense molecular clouds , and useful as a m easure 
of the ionization degree of the gas ( Phillips et al.lll992l) . Since 
H30^ formation requires H2O to exist in the gas-phase, the 
H30^ molecule acts as a natural filter to select hot {T > 100 K) 
molecular gas, and therefore traces more specific regions than 
molecules usually surveyed towards other galaxies. The chem- 
istry of this filtering is the evaporati on of icy grain mantles at 
r ^ 100 K (Ivan der Tak et al.ll2006bh . or that H2O forms in the 
gas phase at high (T > 300 K) temperatures. 

Observations of H^O^ emission will therefore help probe 
the location of dense, warm and active gas in galactic nuclei. 
Combined with information on the H2O abundance, these ob- 
servations may further be used to trace the ionization rate by 
cosmic rays (produced in supern ovae, i.e., starburs ts) and/or X- 
rays (from an AGN). Studies bv Ivan der Tak et aP ([2006a ) have 
demonstrated this use of H30^ for Sgr B2 in the Galactic Centre. 
Furthermore, recent Herschel observations of Galactic sources 
show 984 GHz Oq - 1q H3 0^ absorption in t he diffuse gas to- 
wards GlO.6-0.4 (W31C) dGerin et al.|[20Tg) and 1.03 - 1.63 
THz H30^ emission from the high mass star forming region 
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Table 1. Sample galaxies'' 



Galaxy 


Q'(J2000) 


(5(J2000) 




D 




hh mm ss 




Lo 


Mpc 


IC 342 


03:46:48.00 


+68:05:46.0 


6 X 10** Lo 


1.8 


NGC 253 


00:47:33.12 


-25:17:17.59 


2.1 X 10'" 


2.6 


NGC 1068'' 


02:42:40.71 


-00:00:47.8 


2x 10" 


14.4 


NGC 4418 


12:26:54.63 


-00:52:39.6 


8 X 10'° 


27.3 


NGC 6240 


16:52:58.89 


+02:24:03.4 


3.5 X 10" 


107 


IRAS 15250 


15:26:59.40 


+35:58:38.0 


1.12 X 10'^ 


244 


Arp 299 


11:28:33.13 


+58:33:58.0 


8 X 10" 


42 



a) Ref erences: IC 342:'Downes et al.'('1992'); NGC 253:IStrickl and et al. 
20041) : NGC 1068: Telesco et al. (1984); NGC 4418: Ridg wav et al. 



1994); NGC 6240: ,Yun & Carilli (2002); IRAS 15250: 1 Veilleux et al 



19991) ; Arp 299: lAalto et al . ( 1997) 



W3IRS5 dBenz et al.ll2010t) . In general H3O+ abundances agree 
well with expectations from PDR models and - together with 
H20^ and OH^ - provide important further insight into gas phase 
oxygen chemistry. 

Recently we have detected - 1'^ HaO^ in the nearby 
starburst M 82 and the ultraluminous galaxy Arp 220 
dvan der Tak et al.l2008l) . Derived column densities, abundances, 
and H30^/H20 ratios indicate ionization rates similar to or even 
exceeding that in the Galactic Center. In M 82 the extended 
evolved starburst (PDR) is a likely source of this ionization rate 
while, for the ULIRG Arp 220, an AGN-origin (XDR) is sug- 
gested. 

In the XDR and PDR models grain-processing is not taken 
into account since the chemistry is taking place in the gas phase. 
However, a formation route for H30^ involving the evapora- 
tion (or removal by shocks) of H2O from grain surfaces need 
also to be considered. For example, evaporating H2O reacting 
with HCO"^ may provide an important source of H30^ (e.g. 
iPhilhps et al.]|1992l:lvan der Tak & van Dishoec k"2000^. 

We have used the James Clerk Maxwell Telescope (JCMT) 
in Hawaii to observe the J-^ - - 364 H30^ line (upper 
level energy E^- 139 K) in seven starburst and active galaxies 
which cover a range of environments. Our goal is to use H30^ 
as a tracer of gas properties in galactic nuclei and to see if H30^ 
can serve as a diagnostic tool to distinguish AGN from starburst 
activity. In sections |2][3]and|4]the sample, observations and their 
results are presented. H30^ line parameters are presented in sec- 
tion l4.1l and in section l4!2] H30^ column densities and fractional 
abundances are calculated. In sections ISTI and l572l H30^ abun- 
dances in the context of X-ray and UV irradiated models are 
discussed and in l5.3l the potential importance of grain chemistry. 
In the last section. |531 we present a brief future outlook. 



2. The sample 

We have selected a sample consisting of seven nearby luminous 
starburst and AGNs - and one distant ULIRG. The galaxies 
are selected from their bright HCN line emission. From our 



previous experience with extragalactic H30^ we knew that the 
line is weaker than the standard high density gas tracers such 
as HCN and HCO^ so we restricted ourselves to a relatively 
small sample of seven objects (coordinates, FIR luminosities 
and distances are presented in TablelTJ: 

IC 342 is a nearby barred Scd galaxy of moderate luminosity 
(central 400 pc has LpiR of 6 x 10^ L©) with a central starburst. 
Within its central 300 pc (30") two mo lecular arms end in 
a clumpy central ring of dense gas (e.g. iDownes et al.l Il99"2h 
which suiTounds a young star cluster. The ring is suggested 
to outline the X2 orbits in a larger-scale bar. The chemistry 
of IC 342 has been investigated in detail in a high-resolution 
study by Meier & Turner (2005). They find that the chemistry 
in the ring is a mixture of PDR-dominated regions and regions 
of younger star-forming clouds. The chemistry in the bar/arms 
is do minated by shocks as shown by CH3OH dMeier & Turneii 
12005 ) and SiO dUsero et al .112006 1). Five (A-E) giant molecular 
clouds (GMCs) are found in the ring and arms. The 13" JCMT 
beam of our H30^ observations is pointed towards the region 
of GMC B - but also includes GMCs A and E. GMCs B and 
C are the sites of young ( a few Myr) star formation and are 
also the regions where incoming molecular clouds meet the 
ring. GMC A has more PDR characteris tics. The dust tem - 
perature of IC 342 is estimated to 42 K (e.g. lBecklin et al.ll98dh . 

NGC 253 is also a nearby barred galaxy located in the 
Sculptor group with a compact nuclear starburst and a IR 
luminosity that appears to originate in regions of intense 
massive s tar formation w ithin its central few hundred par- 
sees ( Strickland et al. 2004). From their 2 mm spectral scan 
Martin et al, (2006) suggest that the chemistry of NGC 253 
shows strong similarities to that of the Galactic Center molec- 
ular region, which is thought to be dominated by low-velocity 
shocks. High resolution SiO observations show bright emission 
resulting from larg e scale shocks as well as gas entrained in 
a nuclear outflow dGarcfa-Burillo et al ] l2000l) . Note also that 
it is suggested that NGC 253 is a galax y in which a strong 
starburst and a weak AGN coexist (e.g. IWeaver et al.l 120021: 
Miiller-Sanchez et al. 201 0). High resolution observations of 
HCN and HCO+ 1-0 dknudsen et all |2007|) show strongly 
centrally concentrated emission. The 13" JCMT beam covers 
the bulk of the nuclear emission from HCN and HCO^. The 
central dust temperature of NGC 253 is estimated to 50 K 
(.Melo et al...2002,) . 



NGC 1068 is the nearest example of a type 2 Seyfert galaxy 
luminous in the infrared. Surrounding the AGN there is a 4" cir- 
cumnuclear molecular ring or -disk (CND) and on a larger scale 
there is a NIR stellar bar 2.3 kpc long. This bar is connected to a 
large-scale, molecular starburst ring that contributes about half 
the bolometric luminosity of the galaxy (e .g. Telesco et al. 1984 



ScoviUe et alJ [T988t Flicconi et all [19941 iHelfer & Blitz .1995 



Tacconi et al Bright HCN 1-0 line emission is observed 



towards the CND while the HCO^_J-0 emission is relatively 
fainter by a factor of 1 .5 (e.g. lKohno et alfc OOl). Only the CND 
and a fraction of the bar is covered by our JCMT beam and we 
adopt a source size of 2" for the CND. The dust of the inner 
4" appears to show a strong temperature gradient - from about 
800 K in the very inner region to 150 - 275 K at a distance of 
0."8 or greater (Tomono et al. 2006). lAlloin et all d2000l) find 
temperatures of about 150 K 200 pc from the nucleus. There is a 
also radio jet from the nucleus which falls into our JCMT beam 
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(e.g. lWilson & UlvestaJl983h . 

NGC 4418 is a luminous, ed ge-on, Sa-type ga laxy with a 
deeply dust-enshrouded nucleus (Sp oon et allioOlh . NGC 4418 
is a FIR-exces s galaxy with a logar ithmic IR-to-radio continuum 
ratio iq) of 3 dRoussel et alJ 120031) . This excess may be caused 
by either a young pre-sup ernova starburst or a buried AGN 
jAalto et al.1 feoOTt [ Roussel et al. 2003; Imanis hi et al., .2004) . 



Unus u ally luminous HC^N line emission dAalto et al.l 12002 
l2007t ICostagliola feAaltol I20T0I) has been interpreted as a 
signature of young starburst activity. Mid -IR intensities are 
indicative of dust temperatures of 85 K (Evans et"^ l2003h 
inside a radius of 50 pc (0."5). The IR luminosity-to-molecular 
gas mass ratio is high for a non-ULIRG galaxy indicating that 
intense, compact activity is hidden behind the dust. Recent 
high resolution imaging of CO 2-1 (Costagliola et al. in prep.) 
indicate a molecular source size of 0."5. 



of the receptors of the array has a beam size of 14" at 345 GHz 
and the sources were placed in the receptor labeled HO^. The 
back end was the Auto-Correlation Spectrometer and Imaging 
System (ACSIS), providing 1.0 GHz bandwidth in 2048 chan- 
nels and a maximum resolution of 0.4 kms In order to opti- 
mize baseline stability, double beam switching at a rate of 1 Hz 
was used with an offset of 120". Weather conditions were op- 
timal, with about 1 mm of precipitable water vapor, leading to 
a typical noise rms of 2 mK at 10 km/s resolution. Telescope 
pointing was checked every hour on the CO line emission of 
nearby AGB stars and always found to be within 2". Data were 
extracted into fits files with the starlink software and analyzed in 
CLASS. Linear baselines were subtracted, line parameters were 
calculated by fitting Gaussian profiles to the spectra. The result- 
ing spectra are shown in Fig.[T] Throughout this paper, velocities 
are in the heliocentric frame and redshift is computed using the 
radio convention. 



NGC 6240 is an infrared luminous merger of two massive 
spiral galaxies. Two AGN/LINER nuclei are separated by 1" 
and the bulk of the molecular gas has piled up between the two 
nuclei dlono et alJi2007l: iTacconi erar.ill999h . Luminous HCO+ 
4-3 emission is also emerging from in between the two nuclei 
where the H2 emission is also located. The mid-IR sources (24.5 
//m) are associated with the two X-ray nuclei - the brightest by 
far being the southern nucleus (r=55-60 K)(Egami et al. 2006). 
It has been suggested that the medium between the two nuclei is 
dominated by starbu rst superwinds from the southern nucleus 
(lOhvama et al. 20001) and t he molecular g as is indeed highly 
turbulent dlonoetdj l2007t ITacconi et al J [19991) . The JCMT 
beam encompasses the entire molecular structure of this galaxy 
and we adopt a source size of 1" for the HsO^ emission. 



IRAS 15250+3609 is a relatively distant ultraluminous 
galaxy - probably a major merger with a dominant bright nu- 
cleus. A tidal feature appears to emerge from the southwest side 
of the nucleus and turns around on the eastern side to create an 
enormous closed ring, 27 kpc in diameter. The optical s pectrum 
is a composite of H II and LINER features ( Veilleu x et alj 
I 999|). IRAS 1525 exhibits deep silicate absorption features 
( Spoon et ani2006h indicating a deeply enshrouded nucleus. Our 
JCMT beam covers the whole optical galaxy. 



Arp 299 is an IR-luminous merging system of two galax- 
ies, IC 694 and NGC 3690. Sti'ong '^CO emission has been de- 
tected from the nuclei of IC 694 and N GC 3690 and from the in- 
terface between t he two galaxies (e.g. lSargent & Scovilldl 19911 
lAalto et al]ll997h . The two nuclei, as well as the western over- 
lap region, curren tly undergo intense star formation activity (e.g. 
iGehrzet al.lfT983h . Bright HCN 1-0 emission is em erging from 
both nuclei as well as from the overlap region (e.g. lAalto et alJ 
[T997I) . Water e mission from NGC 3690 suggests that it also har- 
bours an AGN (T archi et al.ll2007h . To cover all three regions a 
map was carried out. 



3. Observations 



The 364.7974 GHz 7k = 3+ - 2" line of H3O+ was observed 
towards the sample galaxies in February and August 2008, with 
the James Clerk Maxwell Telescope (JCMlQ) on Mauna Kea, 
Hawaii. We used the 16 elements heterodyne array HARR Each 



4. Results 

4.1. Line intensities and line widtlis 

Integrated line intensities, line widths and fitted velocities (he- 
liocentric) for the seven observed galaxies can be found in 
Table |2] p-HaO^ line emission was detected towards IC 342, 
NGC 4418, NGC 253, NGC 1068 and NGC 6240, while Arp 299 
and IRAS 15250 were not detected. 

In general it is found that line widths and shapes (Fig. [1]) 
agree well with those found for HCN and HCO^ in the cen- 
tral r e gions of the systeni s (iNguven et al.l 119921 : lAalto et aP 
120071: Knudsen et al.l 12007|) . The most noteworthy deviation 
from this is NGC 6240 where the linewidth of HCN and HCO+ 
(iGreve et al .120091) is more than a factor of two greater than what 
we find for H30^. A possible explanation for this could be the 
weak H^O^ signal resulting in a moderate signal-to-noise pro- 
file. For NGC 1068 the line width agrees well with that of HCO+ 
4-3 but is some what narrower than what is o bserved for HCN 4- 
3 (202 km s ') (iPerez-Beaupuits et al.ll2009l) . 

For IC 342 an unidentified (U) line was detected in February 
2008 - but it did not appear again in August the same year The 
line is somewhat narrower than the H^O^ line and appear at a ve- 
locity of 263 km s"' (rest frequency 364.477 GHz). It is possible 
that it is an atmospheric O3 line. 



4.2. Column densities, excitation and abundances 

Since we have no direct information on the excitation of the 
H30^ molecule we make some assumptions that will have to 
be tested in future observations. The critica l density of the H30 ^ 
transition is high (about n^nt = 10 cm" (iPhillins et al.lll992h ') 
and the low reduced mass of the H30^ molecule makes its exci- 
tation very sensitive to radiative pumping by dust. We therefore 
assume that the coupling between the colour temperature of the 
IR emission and the excitation temperature of H30^ holds as 
long a s the dust colour temperature is above 30 K (see discus- 
sion in Ivan der Tak et al . (2008)). 

For each galaxy we ran a RADEX model ( van der Tak et al.l 

l2007i) with temperatures ranging from r=10-500 K. In these cal- 
culations, the volume density is set to an arbitrary large value, so 



' The JCMT is operated by the Joint Astronomy Center on behalf of 
the Science and Technology Facilities Council of the United Kingdom, 



the Netherlands Organization for Scientific Research, and the National 
Research Council of Canada. 
- For more information see: www . j ach . hawaii . edu/JCMT/ spectral_line 
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»«i IC342 




kni/s 200 



0.02 



0.05 




NGC6240 




Galaxy 


In 


Vline 


6V 


J,. 




[K km/s] 


[km/s] 


[km/s] 


[mK] 


IC342 










08-2008:'' 










H3O+ 


0.30 ± 0.07 


39± 10 


75 ± 18 


3.7 


02-2008:'' 










H3O+ 


0.22 ± 0.06 


40 ±5 


39 ± 13 


5.4 


(U-line 


0.38 ± 0.05 


263 ± 2 


30 ± 5 


11.6) 


NGC 253 


9.35 ± 0.18 


221 ± 2 


180 ± 4 


48.8 


NGC 1068' 


0.52 ± 0.05 


1120 ±8 


160 ± 20 


3.0 


NGC 4418 


1.02 + 0.06 


2129 ±5 


171 ± 11 


5.6 


NGC 6240 


0.32 + 0.05 


7179 ± 15 


177 ± 27 


1.7 


IRAS 15250 


< 0.012^^ 






<1.0 


Arp 299 


< 0.020'' 






<2.0 



a) Integrated line intensities and peak line intensities are given in T^, 
Tlrmmb IS 0.65-0.7 and we have adopted 0.7 here for conservative val- 
ues of A'(H30+). b) Data reduction of H3O+ and U-line taken at ICMT 
02-2008 and 08-2008. c) Icr upper limit for a FWHM linewidth of 500 
km/s. d) 1(T upper limit for a FWHM linewidth of 300 km/s. The point- 
ing center of HARP was in between the two nuclei. We added the spec- 
tra in the three emission regions of Arp 299 covered by HARP to pro- 
duce a spectrum with a Icr upper limit of 2 mK. 



H30'^ abundances relative to H2 {XiYi^O^)) are listed in ta- 
blets] We have used CO 2-1 spectra from the literature and used 
a Galactic conversion factor from CO luminosity to H2 mass 
(2.5 X 10^" cm"^ K"' km"' s) to estimate /V(H2). Note that for 
Galactic nuclei (and starbursts) it is argued that this conversion 
factor o verestimates the H? column density by factors of 5-10 
(see e.g. iMartm et al.l (l2010h ). Since we wish to reduce the risk 
of overestimating the HaO"^ abundances we have however still 
adopted the standard conversion factor. For IC 342 and NGC 253 
we have estimated and average A^(H2) for the 13" beam size. For 
NGC 1068, NGC 4418 and NGC 6240 we have used the same 
adopted source size as for HaO^. 



Fig. 1. Observed spectra of HaO^. The intensity scale is in T^, 
not corrected for the JCMT beam efficiency, that at this frequen- 
cies is 0.7. The parameters of the detected lines were obtained 
by Gaussian fits, shown on the plots as dashed lines. 



that the excitation of HjO^ is thermalized. The adopted exci- 
tation temperature is essentially the radiation temperature and a 
radiative excitation with one excitation temperature is simulated. 
We therefore do not use the full non-LTE capacity of RADEX. 

The results are shown in Fig. |2] Each plot has temperature 
(kinetic or radiation) on the x axis and HsO^ column on the y 
axis. The contour line corresponds to the observed brightness 
temperature. This depends on our assumption for the source size. 



4.2.1. Errors in H3O+ column density and abundances. 

Our assumptions on the excitation and spatial extent of of H30^ 
introduce errors in the column density estimates. For IC 342 in 
particular the source size error can be significant. We assume that 
the source fills the beam which gives us an average column den- 
sity for GMCs A, B and E, but the real column density towards 
a particular cloud should be higher For NGC 253 the double 
peaked line shape suggests that the H^O^ emission is emerging 
from both dense-gas peaks in the center and is filling the JCMT 
beam. For NGC 1068 the H3O+ source size in the CND could 
be anything from <1" to 4". We have assumed a source size 
of 2". For NGC 4418 and NGC 6240 the source sizes ai-e ; 1" 
which means that the beam dilution is significant, but no HjO^ 
emission is missed. 
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Table 3. Parameters" for the RADEX Modeling of the H3O 
Line. 



Galaxy 


Source Size* 




MHsO+f fit 


X(H3 0+)'' 




["] 


[K] 


[cm"-] 




IC 342 


13 


40 


1.25 X 10"* 


4 X 10-' 


NGC 253 


13 


50 


2.25 X 10'^ 


1.5 X lO^** 


NGC 1068 


2 


150 


1 X 10'' 


1 X 10^** 


NGC 4418 


0.5 


80 


2.5 X 10'* 


3 X lO^** 


NGC 6240 


1 


50 


5 X 10" 


2.5 X 10^** 



a) The temperature range investigated for all galaxies was 10-500 K, 
but solutions were fixed where T^^ = Tir. The ^(HsO^) range searched 
was 10'- - 10'^ cm"^. b) For a discussion of adopted source sizes see 
section |2] c) The total beam-averaged column density of H3O* is 2-3 
times higher than that of P-H3O* in Fig.|2l because the ortho to para 
ratio of H3O* drops from the high- temperature limit of o/p = I for T > 
100 K to o/p=2 when T < 50 K. Here we have corrected column den- 
sities for an o/p ratio of 1.5. d) IC 342: yCHa) i s estimated to 3 x 10"^ 
cm-2 (from the CO data of lEckart et all ( fl^Oh ). NGC 253: jVCHz) is 
estimated to 1.3 x 10"' cm"^ (using CO data from' Mauersberger et alJ 
( [1993 ). NGC 1068: From Tacconi et al. ( 1994) the A'CHj) for the CND 
is estimated to be 2 - 10 x 10^^ cm"-. If we take the CO 2-1 data 
from Planesas et alj (1989) and assume a source size of 2" we obtain 
A'(H2)=1.35 X 10^3 cm-^ NGC 4418: Column density for Hj for a 0."5 
CO 2-1 so urce size is A^(H2)= 7.4x 10-' cm"^ (Costagliola et al in prep.). 
NGC 6240; Ilonoet"and2007h estimate NQii) to 1 x 1 0^' cm^^ from their 
CO 3-2 data. 



In Fig.|2l the impact of the assumption of the excitation tem- 
perature on the resulting H^O^ column density is illustrated. 
Note the strong dependence of A^(H30^) for Tex ; 60 K. For 
NGC 253 (for example) A^(p-H30+) changes from 9x10''' cm^^ 
when rex=50 K to 7 x 10'^ for rex=25 K. When Tex exceeds 
60 K, however, the temperature dependence on A^(H30^) is al- 
most gone. Thus for all galaxies, except IC 342, we are not likely 
to overestimate A^(H30^) by more than a factor of 2 since the 
adopted T^xl 50 K. It is possible that H30^ could be collision- 
ally excited. In this case the excitation temperature is likely sig- 
nificantly lower than the ones we assume here (beacuse of the 
high critical density). In Fig. |2]it is evident that the A^(H30^) 
would go up considerably with decreasing temperature. With the 
assumed Tex the resulting A^(H30^) is already qute large for all 
galaxies - even larger H30^ columns would be an interesting re- 
sult indeed, but very difficult to explain with current models. The 
excitation can be constrained through observing multiple transi- 
tions and Requena-Torres et al (in prep.) are currently studying 
the 307 and 364 GHz line of H30^ in a sample of starburst galax- 
ies. 

We conclude that the source size uncertainties dominate the 
errors in the column density calculations. To improve future 
A^(H30^) calculations it is necessary to determine source sizes 
through high resolution observations and to obtain information 
on rex(H30^) through multi-transition observations. 

In addition it should be noted that the estimates of the H30"'' 
relative abundances are dependent on the reliability of the con- 
version factor from CO luminosity to H2 mass. 



5. Discussion 

5. 1. H^O^ abundances and PDR/XDR models 

We compare the column densities of H30^ and H2 in Table |3] 
to chemical models of clouds irradiated by either far-UV, photon 
dominated regions (PDRs) or X-ray photons (XDRs). In general, 
we expect the XDRs to dominate in molecular ISMs surround- 
ing an A GN. The XDR and PD R models we use are pre sented in 
Meiierin k~& SpaansI ( l2005h and lMeiierink et al.1 (l2007h . In gen- 
eral, the thermal and chemical structure of XDRs and PDRs are 
different. Larger parts of an XDR can be maintained at high tem- 
peraure and the ionization fraction of an XDR can be up to two 
orders of magnitude higher (x^ = 10"^ - 10 ') than in a PDR. 
In an XDR H2O may form in the gas phase at high temperatures 
(200-300 K) and then react with either H+ or HCO+ to form 
H3O+. 

Source averaged HCO^ abundances are typically ^ 10"^ 
for all detected galaxies excep t for IC 342 . A general result 
of the models presented in van d er Tak et al.l ([2008) is that the 
X(Hi,0^) does not exceed 3 x 10"' in PDRs - even in the 
models with an extremely high cosmic -ray ionization rate of 
^ = 5.0 X 10 s '. H30^ abundances approaching 10 and 
beyond are more likely to occur in X-ray dominated regions and 
we conclude that the XDR scenario fits the X(Ht,0^) values of 
NGC 253, NGC 1068, NGC 4418 and NGC 6240. For IC 342, 
the results are consistent with a PDR, but a model with a high ( 
0/5.0 X 10-'^ s-i. 

5.1 .1 . Water abundances 

The H30^/H20 ratio is an even better probe of the ionization 
rate, and potential nature of the emitting source. It is straight for- 
ward to obtain H30^/H20 abundance ratios as large as 10"^ in 
XDR models, while for PDR models ratios are generally 10"^ or 
less. Herschel has already measured water abundances in M 82 
dWeiB et al] |2pi0) (see also section I5.4l i and results for other 
galaxies will follow soon. 

5.2. Relative H^O^ and HCO^ abundances and XDR/PDR 
models 

iKohno et al.1 (l200lh and llmanishi et al.l(l2004 find HCN/HCO+ 
1-0 line intensity ratios greater than unity in several 
Seyfert nuclei - w here also the HCN/CO 1-- line ratio is 
high. Furthermore, iGracia-Carpio et aP (l2006t ) find elevated 
HCN/HCO+ 1-0 line ratios in ULIRGs. This is often interpreted 
as a sign of an underabundance of HCO^ compared to HCN 
due to X-ray chemistry. (Although it is important to remember 
that an abundance difference cannot be unambiguously deduced 
from a single-transition line ratio.) Underabundant H CO^ in 
XDRs has been proposed in theoretical work by Malone v et al. 
(1996). However, more recent models by Meiierink& Spaans 
(2005) suggest that HCO^ may be underabundant in mod- 
erate column density (A^h < 10^^^ cm"^) XDRs - but for 
larger columns the reverse is true and X(HCO^) generally ex- 
ceeds X(HCN). Here we compare our derived A^(H30^) with 
A^(HCO^) for the sample galaxies to see if large H^O^ abun- 
dances are paked with a particularly low HCO^ abundance, and 
to compare A^(H30^)/A^(HC0^) with the expectations from cur- 
rent models. 

In the XDR/PDR models, the formation of both H3O+ as 
well as HCO^ is mainly driven by reactions with H2, Hj, 
and H3, and destruction by electrons, unless the electron frac- 
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Table 4. HCO+ZHsO+Ratios. 



G3.13.xy 


iV ^rH^w ) 




IC 342 


2-6x 10"* 


0.2-0.6 


NGC 253 


5 X 10'* 


4 


NGC 1068 


4x 10'3 


24 


NGC 4418 


4.7 X 10'5 


5 


NGC 6240 


3 - 4 X 10'-^ 


0.7 



a) IC 342: Estimated from iNguven et aP (ll992h ('Note that a low 
HCO+/H'^CO+ 1-0 line ra tio of 7 suggests that HCO+ columns may 
be higher. NGC 253: From lMartm"e t al. (2006) (6.5 x 10'^ for T^^=12 
K and for 35" source size. Rescaled A'(HCO+) for the 10" by 14" H3O+ 
source size of ( Knudsen et al. 2007)). NGC 1068: We used the data 
provided in P erez-Beaupuits et al.l |2009) and made a rotational dia- 
gram (assuming a 2" source size) finding rex (HCO+)=36 K. NG C 4418: 
Column density for HCO* (based on data in lAalto et al.l ( l2007h and un- 
published JCMT HCO+ 4-3 data) for a 0."5 source size. NGC 6240: 
We esti mate the HCO^ co lumn density from the excitation information 
given in llono et alj ( 120071) and for a source size of 1". 



tional abundance is very low XfllO^'' - 10^^. However, the 
interpretation of an A^(H30^)/A^(HC0^) abundance ratio in an 
XDR/PDR scenario is not entirely straightforward, but the mod- 
els provide some useful limits. For example, in the PDR mod- 
els HsO^/HCO^ ratios are unlikely to become larger than 3 
dMeiie rin k et a l.| 12010). while for the XDR models it is quite 
straightforward to obtain A^(H30+)> A^(HCO+), but the obtained 
ratio is very column density dependent. When the column den- 
sities are small (e.g A^(H) ^ 10^^ cm"^) ratios as large 20 can 
be obtained, but abundances (and associated brightness temper- 
atures) are small, and in order to obtain significant HCO^ and 
H30''" column densities larger clouds are needed. However, when 
increasing the column, the HsO^/HCO^ column density ratio 
slowly decreases, but moderate ratios around 4-5 are easy to 
reproduce with the XDR model. 

5.2.1 . A?(H30+)W(HC0+) in the sample galaxies 

In tableglwe list A^(HCO+) and A?(H30+)/A^(HC0+) for the sam- 
ple galaxies. The largest A^(H30^)/A^(HC0^) value is found in 
NGC 1068 which has a Seyfert nucleus and its inn er few hun- 
dred pc has been suggested to be an XDR (e.g. Usero et all 
12004'; Pei-ez-BeauDuits et alJ 12009; Garci'a-Burillo et al."2010'). 
The large abundance ratio of ^24 suggested for NGC 1068 is 
consistent with (a low column density) XDR. 

For NGC 4418 and NGC 253 XDR models can explain the 
H30^/HC0^ ratios (but we note that the values with errors are 
also within the range of PDR models). For NGC 253 there is 
evidence that both an AGN and a young starburst is present 
while for NGC 4418 the nuclear activity is so obscured that the 
nature of the activity cannot be discerned (see section |2]i. The 
H3O+/HCO+ ratio for IC 342 is consistent with a PDR model as 
is the relative H30^ abundance. 

For NGC 6240 the H3O+/HCO+ ratio is consistent with 
both XDR and PDR models while the H30^ relative abundance 
favours an XDR. It is interesting to note that even if the molec- 
ular gas of NGC 6240 is collected in-between the two nuclei, 
their radiation could still impact the gas. The two AGNs are sep- 
arated by 1".5 which means that the molecular gas is irradiated 



by intense X-ray emission from two directions on only 370 pc 
distance. This scenario could possibly explain both the relatively 
high H30^ and HCO^ abundances. 

5.3. Grain chemistry and fhe formation of H^Cy. 

Above we have found that the XDR/PDR models can potentially 
explain the H30^ and HCO^ column densities and abundances 
we found for the observed galaxies. However, for NGC 4418 
the problem with an XDR explanation is the large column den- 
sities of gas and dust observed towards its c enter (see footnote 
of Tab. where ICostadiola & Aaltol (l20Toh find global HC3N 
abundances similar to those found towards hot cores in Sgr B2. 
Intense vibrational line emission suggest that the HC3N indeed 
exists in warm-to-hot environments in the very center of the 
galaxy. The HC3N abundances are not consistent with either 
XDRs or PDRs. 

Thus, we have searched for alternative explanations to the 
H30^ line emission outside of the XDR and PDR models. 
Elevated HC3N line emission is consistent with the conditions 
in warm, dense shielded gas associated with embedded star for- 
mation. Is it then possible to obtain relative H30^ abundances of 
10"^ without the X-ray chemistry? 

In a scenario where grains are important, the H2O can evap- 
orate o ff the grains at temperatures >100 K and in lPhillips et al.l 
('1992') the destruction of HCO+ through the reaction of H2O 
is found to be important formation process for H30^ when the 
gas is dense and warm. For a water abundance X iH^O) of 10"^ 
and a density n(H2) of lO^* cm^^ iPhillips et all d 1 9921) find an 
H3O+ abundance, X(H30+), of ^ 10'^ and X(HCO+) a fac- 
tor 2-3 lower. It is unclear however, if the A^(H30^)/A^(H20) 
can be greater than 10'^ w ithin the context of their model. 
Ivan der Tak & van Dishoeckl (120001) find the same effect in their 
study of the of Galactic protostar GL 2136. In their Fig. 1 the 
temperature and density structure is presented including calcu- 
lated concentrations of HCO^ both with and without destruction 
by H2O. 

We suggest that for NGC 4418 it is possible that the H^O^ 
abundance and the N(H20^)/N(HC0^) ratio can be explained 
by a model where the gas is warm and dense and water is com- 
ing off the grains to react with HCO^ to form H^O^. This pro- 
cess could contribute to a reduction of the HCO^ abundance in 
the gas phase which may result in a lower HCO^ line intensity in 
any transition. Note that shocks may be responsible for removing 
the water from the grains m aking it available for further reactions 
including H30^ formation. iFIower & Pineau Des ForitsI (l20I0h 
presents models where water abundances are related to the type 
and strength of shocks and there are a multitide of Galactic ob- 
servational studies showing enh ancement of water emission in 
shocks (Nisinietal. 2010; Lefloch et al.ll2010l: Imiiipfler et al.1 
2010; Melnick et al. 200^. Apart from causing H2O to come off 
the grains, shocks can result in efficient gas-phase formation of 
H2O from OH in the shocked high-temperature regime. 



5.4. H2O+, H3O+, andH20 



Wei B et all (1201 Ol) find a (potentially) surprisingly low 
A^(H20)/A^(H20^) ratio of only a few which they attribute 
to H2O evaporating off the grains through shocks, the H2O 
then becomes photodissociated into O and OH. Ion-molecule 
reactions then make H20^ while H2O abundances remain low. 
It is furthermore iiotewo r thy that X(H20^) > Z(H30+) in M 82 
(Ivan der Tak et al.1 120081: IWeiB et al.l l2010h . This is consistent 
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with the low density (n = 10^ cm ■^ ), high cosmic ray rat e 
(> 5 X 10"^^ s-^) model (Model 2) of 'Meiier ink etaP (l2010h . 
This model does howeve r not include any grain processing. 

In lvan der TakI (l2010h a photoevaporation-ionization hypoth- 
esis for the H30^/H20'^/H20 ratio is discussed for M 82 noting 
that testing this theory requires calculation of the photodissoci- 
ation cross-section of HsO^. Whether this hypothesis holds for 
other galaxies, however, remains to be seen: the H20/H20^ ra- 
tio in M 82 may well be unusually low. In the ULIRG Mrk 231 
A^(H20)» A^(H20^) and the H2O+ hne emission feature is c on- 
sistent with an XDR interpretation jvan der Werf et a i]l2oToh . 

5.5. Future studies 

With this study we have confirmed that the 364 GHz line of p- 
H30^ is feasible to detect in starburst and active galaxies and 
that HsO^ fractional abundances are substantial. The HIFI het- 
erodyne spectrometer onboard ESA's Herschel space observa- 
tory is currently obtaining H2O and H30^ data around 1 THz and 
more accurate estimates of the H20/H30^ abundance ratios can 
be determined for many sources improving the understanding of 
the ionization rates of nuclear molecular regions. Determining 
the cource size and excitation of HsO^ is essential to further im- 
prove the understanding of the impact of the nuclear activity on 
its surrounding interstellar medium and to facilitate a more accu- 
rate comparison with models. In the near future, grain chemistry 
and impact of shocks should be added to the interpretation of the 
data. In particular the results for NGC 4418 - and to some degree 
also NGC 253 - emphasizes this conclusion. 
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Fig. 2. On the x axis we find the temperature of the exciting background radiation and on the y axis the column density of H30^. 
For each galaxy, the brightness temperature of the HsO"^ transition is shown in a gray scale map. The thick black contour marks the 
observed value, with one sigma errors are marked by dashed lines. Background temperatures estimated by IR observations are also 
reported and used to derive a column density consistent with the observed brightness temperature. Note that the colunm density fits 
include only the para column. 



